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Abstract 

We stndy the matrix models calcnlating the sphere partition fnnctions of 3d 
gange theories with Af = 4 snpersymmetry and a qniver strncture of a l) 
Dynkin diagram (where each node is a nnitary gange gronp). As in the case of 
necklace {A) qnivers, we can map the problem to that of free fermion qnantnm 
mechanics whose complicated Hamiltonian we find explicitly. Many of these 
theories are conjectnred to be dnal nnder mirror symmetry to certain nnitary 
linear qnivers with extra Sp nodes or antisymmetric hypermnltiplets. We show 
that the free fermion formnlations of snch mirror pairs are related by a linear 
symplectic transformation. 

We then stndy the large N expansion of the partition fnnction, which as 
in the case of the A-qnivers is given to all orders in 1/A^ by an Airy fnnction. 
We simplify the algorithm to calcnlate the nnmerical coefficients appearing in 
the Airy fnnction and evalnate them for a wide class of l)-qniver theories. 
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1 Introduction 


Three-dimensional gauge theories with A/” = 4 supersymmetry have an SU{2)c x SU{2)h 
/^-symmetry group, with SU{2)c acting on the Coulomb branch moduli and SU{2)h 
acting on the Higgs branch moduli. A large class of three-dimensional AA = 4 gauge 
theories flow to interacting fixed points in the infrared limit. These theories possess a 
remarkable duality known as mirror symmetry [1], which is the statement that pairs of UV 
gauge theories flow to the same infrared fixed point with the action of SU{2)c and SU (2) h 
exchanged. In particular the classical (protected) Higgs branch of one theory matches the 
quantum-corrected Coulomb branch of the other theory m- The exact results obtained 
from the technique of supersymmetric localization provide a privileged testing ground for 
mirror symmetry. In particular the full partition function of 3d AA = 4 gauge theories 
defined on a three-sphere reduces to a relatively simple matrix model |1]. This matrix 
model is independent of running coupling constants and therefore computes the three- 
sphere partition function of the infrared fixed point. Mirror symmetry was successfully 
tested by matching the matrix models of pairs of mirror-dual theories in ISHU- More tests 
of mirror symmetry by matching the Coulomb branch and Higgs branch Hilbert series of 
dual theories were achieved in iHII]. 

The most studied J\f = 4 gauge theories subject to mirror symmetry are infrared 
fixed points of quiver theories of type A, D or E, referring to the shape of the quiver as 
an extended Dynkin diagram and likewise to the orbifold singularity in their M-theory 
realization [laiia]. It was found in na that the mapping between matrix models of 
mirror-dual theories of type A can be expressed as a very simple canonical transformation, 
exchanging position and momentum, in the Id free fermion formalism developed in jl5j . 

The free fermion formalism arises from the observation that the matrix model com¬ 
puting the three-sphere partition function of A/" ^ 3 Chern-Simons-Matter quiver theories 
can be re-expressed as the partition function of a gas of non-interacting fermions in one 
dimension with a non-trivial Hamiltonian. This formalism allowed the use of powerful 
techniques from quantum and statistical mechanics to solve the matrix models as an h 
expansion, h being related to the Chern-Simons levels of the 3d theory. It was found 
in [15] that the perturbative part of the partition function at large N takes the form of 
an Airy function. 

This Airy function behavior was first found for the particular case of ABJM theory 
[TB] . in [HI, based on its original large N perturbative solution [TB]. Non perturbative 
corrections have been studied intensively for ABJM [TUII27] . ABJ [2S1I2H], and some other 
examples of A-quiver theories janiisi]. Recently, the grand partition function of ABJ(M) 
was determined exactly |32l|33], in perfect agreement with numerical results 131. 

In this paper we find the free fermion formalism associated to infrared fixed points of 
A/" = 4 Yang-Mills quiver theories of type D and similar Chern-Simons-matter theoriesjl] 
A Z)L+ 2 -quiver is described by a linear chain of (L — 1) U{2N) gauge nodes, connected by 
bifundamental hypermultiplets, with pairs of U{N) nodes attached at both ends of the 
U{2N) chain, as shown in figure [U In addition each node may couple to an arbitrary num¬ 
ber of fundamental hypermultiplets, indicated by the boxes in the diagram. Contrarily 
to type A theories, mirror symmetry does not relate pairs of type D theories. Instead, 
i)-quivers are mapped under mirror symmetry to linear quivers with U(2N) gauge nodes 
and extra symplectic gauge nodes Sp{2N) or antisymmetric hypermultiplets at both ends 

^Chern-Simons terms generically lead to only Af = 3 supersymmetry, and we mostly set them to zero 
after section [21 
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Figure 1: The general -DL+ 2 -quiver with arbitrary fundamental matter. 


of the quiver chain [T^[35H57] . that we call generically linear quivers^ 

Roughly speaking, the length of the linear quiver mirror depends on the number of 
fundamental hypermultiplets coupling to the U{2N) nodes and the structure of the end 
of the linear quiver depends on the number of fundamental hypermultiplets on the U{N) 
nodes. When the number of fundamental hypermultiplets on the two U{N) nodes on one 
side of the Z)-quiver are equal, the mirror theory has a terminating Sp{2N) node. When 
the numbers of fundamental hypermultiplets differ by one, the mirror theory has a termi¬ 
nating antisymmetric hypermultiplet. For instance if the numbers of U{N) fundamental 
hypermultiplets are given by = 0 and = 1 (in the notations of 

hgure d]), then the mirror theory has Sp(2N) nodes at each end as in hgure [2al while 
nio) = ^'{0) = ri^iL) ^ g ^/{L) ^ ^ leads to a mirror linear quiver with a terminating 
Sp(2N) node at one end and an antisymmetric hypermultiplet at the other end as in 

hgure ESI 

When the numbers of fundamental hypermultipets on the two U{N) nodes differ by 
two or more, we do not know what is the precise mirror dual theory. Naive considerations 
on the Hanany-Witten IIB brane setup realizing these quivers [3611^ . lead to ‘bad’ mirror 
dual linear quivers m- ‘Bad’ quivers are quivers whose gauge group cannot be completely 
higgsed and whose three-sphere partition function is divergent. The formal manipulations 
leading to a Fermi gas description presented in section El can be done for any quiver, 
however when we study mirror symmetry, we focus only on duality between ‘good’ quivers. 




(b) 



Figure 2: Examples of linear quivers that are mirror dual to Z)-quivers 

^There are many other types of linear quivers, which we do not study in this paper. We use this name 
to refer to the quivers as in figure [2] and hope that this will not cause confusion. 

^Note that the parameters nF) in figures ITI and [2al I2bl do not get mapped to each other under mirror 
symmetry. The actual mirror map is more involved, see section [S] 
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Our starting point is the matrix model computing the three-sphere partition function 
of the quiver theories. This matrix model is expressed as an integral over the Cartan 
subalgebra of the gauge group and the integrand is a product of classical contributions 
and one-loop contributions in the (exact) saddle point analysis of supersymmetric local¬ 
ization [1] 




|W| 


dA ^r'lfl.RS * * Z] 


hyp 5 


Cartan 


( 1 . 1 ) 


where |W| is the order of the Weyl group. The classical contribution depends on Fayet- 
Iliopoulos (FI) and Chern-Simons (CS) parameters. These are given by 


7 FI ^ 27riCTrA 7 CS ^ vrifeTrA^ 

■^class ^ ) -^class ^ ) 


( 1 . 2 ) 


with C the FI parameter, k the Chern-Simons level and Tr the trace in the fundamental 
representation. The one-loop contributions of the AA = 4 vector multiplet and hypermul- 
tiplet in a representation TZ are 

^vec = n ^ sinh^(7r a • A), Z^yp = ^ , (1.3) 

Lo l^^2cosh{n{wX + m)) 

where a runs over the positive roots of the Lie algebra and w over the weights of TZ. m 
is a real mass parameter for the hypermultiplet. We provide in appendix the explicit 
matrix factors relevant to the H-quivers and their mirror linear quivers. 

By manipulating the matrix models of H-quivers we are able to re-express it in the 
form 



where p depends on the content of the theory and Ua is the number of cycles in the 
permutation a. This expression differs from the analogue expression for A quivers by the 
presence of the factor 1/2”'^. Dehning the density operator p by (A| p |A') = p(A, A'), we 
are able to recast the partition function of theories with vanishing mass and FI parameters 
into 

^*=42 (-1)” f ■ (1-5) 

■ (ts5jv j=l *=1 ^ ^ 

where R is the reflection operator i? |A) = |—A). This can be interpreted as the partition 
function of N non-interacting fermions on a half line with Neumann boundary condition 
at the origin, with a Hamiltonian H = — log p. We also hnd that can alternatively 
be recast into the partition function of N non-interacting fermions on a half line with 
Dirichlet boundary condition at the origin, which translates into having the projection 
instead of in fll.5p . This freedom exists because, as we show in appendix O 
the density operators associated with D quivers have the remarkable property that their 
spectrum is pairwise degenerate between even and odd states on the line. 

Similarly we hnd that the matrix models associated to linear quivers can be set in the 
same form 01.411 or 01.511 . with the same property that the density operator has pairwise 
degenerate spectrum. The derivation of the density operators associated to the D-quivers 
and their mirror-dual linear quivers are presented in section [21 

Having found the density operators for £)-quivers and linear quivers, we observe in 
section [3] that density operators of mirror theories are the mapped by the transformation 
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on position and momentum operators 


q^-p, ( 1 - 6 ) 

as was first found in [H] for type A quivers 0 

With a free fermion formalism, we are in a position to evaluate the perturbative part of 
the large N expansion of the partition function of -D-quivers, as was done for A-quiver 
theories in [39]. We present the analysis in section 0] and derive explicit expressions for 
theories with vanishing masses and FI parameters. This is conveniently described from 
the grand potential which is the logarithm of the grand canonical partition function 


00 

S(z) = 1+2] Z{N)z^ = z = e^^. 

N=1 

Using phase space techniques, the grand potential evaluates to 

C 

J{fi) = —pA + Bfi + A + 0{e~°‘A ) a > 0 , 

o 

which after inverting the relation fll.7p leads to 


Z(N) 


C-ie''Ai 


C-i{N-B) +ZUN), 


(1.7) 


( 1 . 8 ) 


(1.9) 


where Ai is the Airy function and Z^ZN) denotes non-pertnrbative, exponentially snp- 
pressed contributions at large N. For quiver theories without CS terms and without mass 
or FI deformations, we are able to hnd the exact values of the coefficients C and B, 
but not A, nor the non-pertnrbative corrections. The method we develop is simplihed 
compared to earlier literature. 

The fact that the result for H-quivers is similar to that of A-quivers is consistent with 
the localization calculation of supergravity in AdS^ [ID] . That suggests a universal answer 
for all conformal 3d theories with enough supersymmetry and an M-theory dual, which 
includes both the A and D-quivers. 

The free fermion formalism for 3d A/" = 4 theories with no U(2N) and a single Sp(2N) 
or SO{2N) node was studied previously in [H]. We compare our formalism to theirs in 
appendix [PlFI 

We conclude with a discussion of some open questions and possible extensions of our 
work in section (5] We also include a simple holographic test of our results. 

During the course of the work we learned that similar questions were studied by 
Moriyama and Nosaka, whose results are published concurrently 021 . 


Notations: We consider quivers with nodes of rank N or 2N and we use the indices i,j 
when the label runs over 1, ■ ■ ■ , N, and /, J when the label runs over 1, ■ ■ ■ , 2N. Moreover 
we dehne 

sli X 

shX = 2sinhTTx , chx = 2 cosh ttx , thx = —— = tanh7rx. (I-IO) 

chx 

^Mirror symmetry for type T-quivers extends to an SL{2,'Z) gronp of dualities [71131]. We are not 
aware of a similar extension for Zt-quivers and consequently our discussion concerns only the mirror 
transformation corresponding to the so-called S'-duality. 

®We point out a mistake in their next-to-leading contribution to the grand potential. 
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2 D-quivers and linear quivers as free Id fermions 

In this section we show how to re-express the partition fnnction of D-qnivers and linear 
quivers as the partition function of a gas of free fermions on a half line. 

2.1 Free fermions formalism for A quivers 



Figure 3: An Ai_i-quiver. denotes the number of fundamental hypermultiplets at the 
corresponding U{N) node. 


We start by reviewing how the free fermion formalism arises for J\f = 3 A quiver 
theories with U{N) gauge groups [I5]. The A^^.i-quiver diagram is presented in figure [3l 
The matrix model for these theories is given generically by 


Z{N) = 




r Yi n 

a^l i 






l(a)^ 


LKJ 


Ui,j ch(Al 


(a) 




(a+l) 


-1- ?7i(“)) 


( 2 , 1 ) 


where i = I,-- - ,iV are the eigenvalues of the node (with the identification 
^(^+1) = and are the bifundamental masses. The factor n. ^(a)(^(a)) represents 
all terms which depend only on single eigenvalues and includes contributions from the CS 
term with level , the FI term with parameter and fundamental hypermultiplets 
with individual masses AA ,« = !,•••, 


9,(“) 


p{a) = g7rjfc(“)A2 g27riC(“)A 


A ch(A + /ia^) 


N 


Using the Cauchy determinant identity 

ni<j sh(Ai — \j) sh(Ai — Xj) ^ (—1')'^ rr 

Ui,j ch(Ai - ~Xj) \A ch(Ai - A^P)) ’ 

the partition function can be re-expressed as a sum over L permutations 


Z{N) = 


N\^ 


2 

ct(“)sSjv “=1 * 


irU)(A(“)) 


ch(Al“> - + m(«)) 


( 2 . 2 ) 


(2,3) 


(2.4) 
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By relabelling the eigenvalues one can factor out all but one of the permutations, picking 
up an overall factor of This gives 


. ^ L N /L-1 

cteSn “=1 * Va=l 




ch(AS“^ - +m(“)) 

F(i)(A<'-)) 


(2.5) 


ch(A<''’ - A A) + ">('-)) 


This integrand is a series of kernels of pairs of specihc eigenvalues of successive nodes 
ultimately coupling each Awith . This can be encoded graphically by the following 
diagram 

(A«} ^ {A(2)} - ... - {A(^)} - {A«} . (2.6) 

One can express the kernels in terms of canonical position and momentum operators 
q, p, which satisfy [g,p] = ih. Taking A to be position eigenvalues, we have 

F(A)5(A'- A) - <A'| F(,)) |A> , = <A| ^ |A'> , |A> = |A - m> . 

(2.7) 

Here we used the standard relation between the position and momentum bases \p) = 
5 have (All A 2 ) = 5(Ai — A 2 ) and (pi| P 2 ) = ^{Pi — P 2 ). We choose to 

normalize p such that ^ This allows one to write the integrand of fl2.5p as 




^27rim^^^p ^27Tim^^^p 

' |Af'><Af>|F«(a.' 


chp 


chp 


Af)>...(Af)|FW(4) 


chj3 


We obtain the hnal expression for Z{N), 



( 2 . 8 ) 


Z(N) 



(2.9) 


with 


P = 




clip 




chp 


F<‘>(«) 


chp 


( 2 . 10 ) 


This expression coincides with the partition function of N non-interacting fermion^ living 
on a line, with a Hamiltonian H given by p = e~^. In this language. A* is the position of 
the fermion on the line. 

Clearly fl2.9p is fully determined by the spectrum of p. In this formulation there is thus 
a natural splitting of the computation of the partition function into two distinct steps. 
The hrst is to hnd a suitable density operator that encodes the desired quiver theory, 
and the second is to solve the resulting quantum mechanics problem. In this and the 
next section we concern ourselves only with hnding the density operators and studying 
relations between them, we deal with the explicit computation of the partition function 
in section 0] 

For the remainder of this section we suppress the hats above operators. 


®The fermionic statistics is understood from the anti-symmetrization over permutations of the positions 


A,. 
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2.2 D-quiver analysis 

We now implement the Fermi-gas formulation to the matrix models computing the 
partition functions of Z)L+ 2 -quiver theories with unitary gauge groups. These theories are 
characterised by gauge group U{2N)^~^ x U{NY and a quiver structure as in hgure [H 
We work out the full details in a simple example, a T) 4 -quiver theory with n fundamental 
hypermultiplets attached to a single U{N) node. We do not turn on mass or FI defor¬ 
mations. The quiver diagram for this theory is shown in figure IH The generalization to 
other D^-quivers is outlined in section 12.41 



Figure 4: A l) 4 -quiver with n fundamental hypermultiplets on one of the U{N) nodes 


Our aim is now to hnd a suitable density operator that repackages the matrix model 
for this theory into an expression like 02.91) . In order to do this it is useful to collect the 
eigenvalues of pairs of terminal U{N) into a single set of 2N eigenvalues. Anticipating 
this, we label the eigenvalues of one pair of U{N) nodes respectively bjEl and 

^2") f2') 

i = 1, ■ ■ ■ , N, and the other pair by A^ y The eigenvalues of the U{2N) node are 

labelled ■, / = 1, • • ■ , 2N. Following the rules from appendix the matrix model for 
this theory is given by 


Z{N) = 


r Yl Yl P? - Pp) 

'J n — n 7 


m*(2N)\ 

n..j (4°’ - Af) (A'* - >■?) - Agt,) 

ni,j d' - -'A) d' (.'w+i - 'll (-'f’ - XX d' (^w+i - X’) Hi cii" 4'^’ 

( 2 . 11 ) 

Including artihcially the factors ni,,sh(A<“'-AW,) n, j sh [xP — XP_^j) in the numerator 
and denominator, one can use the Cauchy identity 02.31) . as well as a modified version 


Yli<j “ P) , , s MAA) 

-^--- = (-1) ^ 


N 


o-g 5 jv *=1 '^o-p)) 


Ylij sh(Ai - Xj) 

to re-expressed Z{N) as 

r7(AT\ 1 'W ( iw(0)+a(2)+r(0)+r(2) frT ,2VWa) 1 

r(“)G52jv 


( 2 , 12 ) 


N 

n 


2N 

n 


d =h (A'"' - A™ =h (AP - A^y„,„) U ch (AP - A';i„„) ch (AP - A«,,„) 

(2.13) 


^This is a slight abuse of our earlier notation where superscripts distinguished gauge group factors. 
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Successive relabelings of the indices allow us to remove the sum over and 

and compensate for it by an overall factor of A^!^(2A^)! 


Z(N) = 


NP 


N 


f, 2 N 

s (-i)d 

rsS2JV «=0 *=1 


h! ch'^ 


2N 


n 


n 


U sh (A<"' - A™ i) sh (a!"> - Ag’^J U ch (xf > - A«) ch (A« - A^,*,) ' 

(2,14) 

As in the case of the A-quivers fl2.10p . we would like to write this as the successive 
interaction between pairs of eigenvalnes. Dehning the reflection permutation R by 


R{i) = N + i, R{N + i) = i, 


(2.15) 


the integrand of the matrix model can be viewed as a series of kernels pairing eigenvalues 
of adjacent nodes in a chain that goes back and forth along the quiver, according to the 
diagram (c./-, 112.6^ 1 

{A(2)}oh. (2.16) 

r ^ 

Traversing the quiver back and forth we end up with the composite permutation 

Rt-^Rt, (2.17) 

so we can write the partition function in terms of a kernel relating A^^^ and 
Note however that another eigenvalue of the central node is integrated over to 

get this kernel. So for each permntation r we need to choose half the eigenvalnes of A*^^^ 
on which the kernel acts. Let us call the set of N indices of those eigenvalues JC{t). It 
is chosen to be closed under the composite permutation Rt~^Rt and such that R takes 
this set to its complement R{K{t)) = JC{t). The partition function can be expressed in 
the following way 


Z(N) = 


NP 


eS2N “=0 i=l keK.(T) 


TeS2N 

_iy{T{k)) 


ch”Af> ", ch(A“-A® ) 


(Art, - Agh)) * (Agh) - KRrm) <=‘‘ (A- At“>.H,(,)) 

Aio) TTw TTw V ’ 

\\-^RT{k) '^RT~^RT{k)) '^Rr-^Rrik)) 


(2.18) 


where 


s{k) 


0, k = !,■■■ ,N , 

1, A: = iV + l,--- ,2iV. 


(2.19) 


To be able to write the partition function in terms of a density operator we need to 
include the contribution from the fundamental hypermultiplets into the product over k 

in fl2.18p . However, the fundamental hypermultiplets couple only to the eigenvalues A'^^ 

(2) (2) 

with i = !,■■■ N, which depending on r is either or but not both. These 

















two options happen with equal probability for each combined permutation Rt so 

we can write it as the sum (normalized by 1/2^)11 


Z{N) = 


2 ^iV !2 

n 

keKir) 


Y, (- 1 )' fn<i“AW n 

tgS2n a—0 /cG/C{r) 

1 


ch sh (A®) - Ag(;-)) ch 

1 


+ 


RT{k) V ^_Rr(fc)> 

1 


ch - XZkd sh (aS(,) - ch (A^) - Ay.^^^y J ch- A^^.^ 
1 1 1 


d 2 ) 


Ch (A^y^ - A 


-iRT(fc)) sh (A^-y^(fc) ^nl-iRTik)) ch (A 


( 0 ) _ Ri) 

fiT-lijT(fc) '^RT~^RT{k) 


— \ ’ 
Ad _-1 

( 2 . 20 ) 

where we have used nA:sA:(T)(“^)^^^ RRk)) = "hhis expression can now 

be recast as 


2(^) = «ivW E (-1)" n (2.21) 

• rsSzjv 


with 


P(A,A') = 2 


Jn 

a=l 


dXr. 


keK(T 


-1 


1 1 
+ 


ch(A - Ai) Vch- Ai sh(Ai - Aa) sh(Ai - Aa) ch- Aa 
fill 
ch(A 2 - As) ch(A3 - A 4 ) sh(A4 - As) ch(A5 - A') ’ 


( 2 . 22 ) 


where we chose the normalization factor for convenience. 

The kernel p dehnes a density operator through the relation p(Ai,Aa) = (Ai|p|Aa), 
which has a representation in terms of canonical position and momentum operators fl2.7p 


P 


1 1 
2 chp 


shp ^ shp 


ch- g ch p ch p ch- q 


shp 
ch'^p ’ 


where in addition to fl2.7p we have used 


1 

sh(A - A') 



shp 

chp 


(2.23) 


(2.24) 


To make further progress, we need to study the combinatorics of the composite per¬ 
mutations Rt~^Rt. We relegate these additional technical calculations to appendix [B] 
and provide the hnal simplihed result, which involves only a sum over permutations of 
Sn 

= iv! E ri"- . ( 2 - 25 ) 

creSff i=l 2=1 

®In more detail, note that Rt~^Rt remains the same if one multiplies r on the left by any combinations 
of two-cycles appearing in R. For a given r this generates a set of 2^ terms, half with T(fc) in {1, • • • , N} 
and half in the compliment. 
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where is the number of cycles in a. 

Because of the factor ^ fl2.25p cannot be interpreted directly as the partition 
function of N non-interacting fermions. When all FI and mass parameters are turned off, 
we hnd that it can be understood as resulting from a projection onto half of the states of 
a fermionic system. We show in appendix [C] that the density operator p commutes with 
the reflection operator i?, dehned by 

^|A>=|-A>. (2.26) 

Consequently, the Hilbert space can be split into even and odd eigenstates. Furthermore, 
the spectra of even and odd eigenstates are identical, allowing us to rewrite Z{N), using 
the projector as 


1 r ^ ^ 

Vi iA,(.)>. (2.27) 

■ o-eSjv j=l j=l ^ ^ 

This can be readily interpreted as the partition function of N non-interacting fermions 
at positions |Aj| on a half-line with a Hamiltonian H = — log p where the operator is 
responsible for the projection onto particle states with even wavefunction on the line, or 
equivalently particle states on a half-line with Neumann boundary conditions. 

Likewise we can use the projector to express the partition function in terms of 
the odd states 


Z(N)- 


N\ 


^ N N 

^(-ir 

o-eSn i=l i=l 


1 - R 


I Acr(i)) 


(2.28) 


In this case we would interpret Z{N) as the partition function of N non-interacting 
fermions on a half-line with Dirichlet boundary condition at the origin. 

We did not hnd a free fermion interpretation of the partition function (12.251) for the 
cases with non-vanishing masses and FI parameters. 


2.3 Linear quiver analysis 

We turn now to studying the matrix models of U{2N) linear quivers terminating at each 
end with either an Sp{2N) node or an antisymmetric hypermultiplet. In many cases these 
quivers are known to be the mirror duals of Z)-quivers [2]. 



Figure 5: A linear quiver with an Sp{2N) node at one end and an antisymmetric hypermultiplet 
at the other. 

We perform the analysis in the case of the quiver with a single U(2N) node which 
is coupled to a single Sp(2N) node and has also an antisymmetric hypermultiplet. This 
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example contains all the ingredients to treat any other qniver of this class, as we elaborate 
on in section 12.41 

The Sp{2N) node in onr example has eigenvalnes labelled i = 1, ■ ■ ■ , N and 
fnndamental hypermnltiplets. The U{2N) node has eigenvalues I = 1, • • • ,2N an 
antisymmetric hypermultiplet as well as fundamental hypermnltiplets. To simplify the 
expressions we do not include any FI terms or masses. Following the rules in appendix |3 
the matrix model is given by 


ZiN) = 


2''Af!(2Af)! J 0 U ,^(0) U ;^(i) 


(Af - Xf>) sh= (A|»> + Xf>) n,<,,sh= (A 


n^cipAf-A'‘’)ch(Ar+A<‘’)n 


,( 0 ) 


(0) 


(0) 


AO)' 

J_ 

di)> 


db 




(2.29) 


i<j 


ch (aW + a“) 


A hrst step is to write the contribution of the Sp node in terms of 2N eigenvalues 
satisfying 

= -^1°’ ■ (2-30) 

The interaction between the Sp(2N) and U{2N) nodes combine to a single Cauchy de¬ 
terminant 02.31) 


(Af - Af) sh^ (Af + Af) nf=iSli2Af' n,<2sh (Ai‘> - A<‘') 


do) , \(0) 


do) 


db x(b' 


n 


i<j 


n„ch(Af>-A<'')ch(Af'+A«) 
sh (Af> - aJ'') sh (A“ - A™) 


2N 


n„ch(A<»>-A3') 


dbA 


2 ^ n ^(1) A 

(0)sS2jv -f=l ™ “ \(0)(I)) 


(2.31) 

The remaining terms involving the eigenvalues of the U{2N) node can be interpreted as 
a Pfaffian, rather than a determinant. We can use the identity pi 143] 


n 


Xj — Xj 

1 -h XjXj 


Pf 

\l + XjXjJ 


1 

2^m 


N 

s (-irn 

tGS2N ^=1 


1 “h ^T{i)'^TR{i) 


(2.32) 


where R is again the permutation R{i) = N + i modulo 2N. Plugging x = we 

obtain 


n 


(aP - x'j>) 

/<j ch [X\ ^ + Aj^) 


N 




(b 


2^A^! 


S (-irn 


TeS2N 


■=i ch 


(b 


(2.33) 


As before, we can remove one of the permutations coming from 02.3ip and 02.33p by 
a relabelling of eigenvalues, picking up an overall factor of (2iV)!. This gives 


Z(N)- 


1 


22W fjl2 


2 (-1)" f d'‘A<»><i“A<‘> n 

C Q.-. » r i 1 

1 sh (A 


^ sh2Af) 


'reS2N 
N 


dl ch2"^°’ Af ^ jj, ch"^^' A 


n 


,n(l) x(b 


db 

r(i) 




) 


(2.34) 


n ch (At»> - At;;,) ch (a;;;, + a™,,) ch (a™ - a™ .,) 
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replacing for convenience t t we can again rewrite the expression as a product over 
the set /C(r), consisting of N indices closed under the permutation Rt~^Rt 


Z{N) = 


22VjY!2 


2 (-1)" [ 

rs52jv keK(i 


(l)dfc) sh2Af 


1 


(_!)•«»)) 


Af ch""’ aP ch"'” A« , 

,( 1 ) 


ch (A™ - A™) 


( 1 )' 


ch (A<‘' + A<‘J. 


RT(k) 


) 


cl> + A 


( 0 ) 

RT{k) j 


(2.35) 

The (—l)*h(^)) signs comes from re-expressing sh in terms of k, while 

fl2.30p is responsible for the (—signs as well as allowing the replacement in the last 
denominator 


^(0) =_\(o) 

T~^ RT{k) Rt~^ RT{k) 


(2.36) 


As in the case of the T)-quivers fl2.2ip . we obtain a density operator between two A^°^ 
eigenvalues related by the permutation Rt~^Rt 


p(A,A') 



sh2A 1 1 sh(Ai — A 2 ) 1 1 

ch2-'°^ A ch(A - Ai) ch"'^' Ai ch(Ai + A^) ch"'^' A 2 ch(A2 + A') ’ 


(2.37) 


in terms of which the partition function is given exactly as in fl2.2ip . Expanding sh(Ai —A 2 ) 
and reversing the sign of A 2 allows us again to represent the operator in terms of canonical 
position and momentum operators 


P 


1 sh2g 1 1 

2ch2"™gchpch"^'^g 


1 sh2g 1 1 

2ch^n(o) qchpclR^^^ q 


(sh q —— ch g -I- ch q —— 
\ ch p ch p 

\ chp chp 



I —--— 

ch'^^'^chp 
1 1 
ch^'^^chp ■ 


(2.38) 


The same arguments as for the D-quiver allow us to express Z{N) as the partition 
function of N non-interacting fermions on a half line with Neumann fl2.27p or Dirichlet 
02.281) boundary conditions. Indeed, as we show in section [21 the density operators of many 
linear quivers are related to those of iA-quivers by a linear canonical transformation. 


2.4 Generalization to longer quivers 

It is straight-forward to generalize the analysis to quivers with an arbitrary number of 
U{2N) nodes and arbitrary number of fundamental hypermultiplets on each node. Just 
as for the A-quiver theories, the matrix model contributions from hypermultiplets trans¬ 
forming in the bifundamental representation of pairs of U{2N) gauge nodes combine with 
the vector multiplet contributions to form one Cauchy determinant between each pair of 
adjacent nodes. This translates into ch“^p terms in the density operator. We represent 
the contribution from fundamental hypermultiplets and the FI and CS terms of the 
node again by F*^“^(g) fl2.2p . This leads to a piece in the density operator of the form 


1 

chp 


F>»{q) 



1 

chp 


F(^\q) 


1 

chp 


(2.39) 


In all of our examples the density operator combines kernels going back and forth along 
the quiver. The contribution from the ends of the quivers {U{N) nodes for Z)-quivers and 
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S'p(2A^)/antisymmetric hypermultiplet for linear quivers) are the same as in the previous 
sections. 

For the D-quivers the contribution from going back along the quiver as in fl2.16p gives, 


J-F<3)(,)-LF<=)(5)-i-FW(,)-L 

ch p chp ch p ch p 


(2.40) 


For the linear quivers, the antisymmetric hypermultiplet or Sp node introduce a minus 
sign, like the replacement A 2 ^ — A 2 that gave fl2.38p from fl2.37p . Therefore the second 
part of the density operator includes 


1 

chp 




1 

chp 




1 

chp 


F^i-q) 


1 

chp 


(2.41) 


2.4.1 Generalized D-quivers 

Let us consider a Dl +2 quiver with arbitrary number of gauge nodes and fundamental 
hypermultiplets on each nodej^ as shown in hgure [H We label the U{2N) nodes by 
1, • • • , L — 1, and as in hgure [T], we distinguish parameters from pairs of terminal U{N) 
nodes by primes: F^^\ p''io)^piL)^ F'^^\ We note that all bifundamental hypermultiplet 
masses can be set to zero by shifting eigenvalues^ The above rules lead to the density 
operator 



(2.42) 

We can easily recover fl2.23p by setting L = 2, F^'^'> (q) = ch“"' q and all other F^^'> = 
F'F) = 1. 


2.4.2 Generalized linear quivers 

We can proceed in a similar fashion to write down the density operators for longer linear 
quivers, where each end of the U{2N) linear chain has either an Sp(2N) node or an 
antisymmetric hypermultiplet and any number of fundamental hypermultiplets on all the 
nodes0 Again we note that the masses for all bifundamental hypermultiplets between 
U{2N) nodes can be set to zero by shifts of the eigenvalues. We cannot always do the 
same for the mass of antisymmetric hypermultiplets, or for those charged under Sp{2N), 
so we keep these masses as well as the fundamental hypermultiplet masses. We consider a 
quiver with L U{2N) nodes and again package the FI, CS and fundamental hypermultiplet 

®Note that in order for the matrix model to be convergent, such a theory must have at least one 
fundamental hypermultiplet. With this condition violated the formal manipulations still go through, but 
the divergence of the matrix model translates to a density operator which is not of trace class. 

^*^Shifting eigenvalues to remove masses can introduce an overall phase in the matrix model. Such phases 
are unphysical, in the sense that they arise from background (mixed) Chern-Simons terms that can be 
added to the regularization scheme when computing the partition function of the 3d theories [441H5] . 

While the formal manipulations again go through in all cases (see also footnote |9]), convergence of the 
matrix model requires a total of at least three fundamental hypermultipets, with at least one coupling to 
the terminating nodes at each end of the quiver. 
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contributions into Two instances (out of three) of such general linear quivers are 

pictured in hgure [2j The density operator is given by 


^L-l 


^L-l 




I ]qF(“)(,)—j FW(,)B(t+«(p,,) I ]q j F<'>(-«), 

(2.43) 

where the functions (p, q) account for whether the ends of the quiver terminate with 
an Sp{2N) node or antisymmetric hypermultiplet. 

If the quiver terminates with an Sp(2N) node we have 


Bsi i.P, q) = 


^27TirrS^^p 


chp 


sh(2g)F(“)(g)- 


chp 


(2.44) 


where is the bifundamental mass and (q) captures the contributions from the CS 
term with level and fundamental hypermultiplets of Sp{2N) with masses 


27rifc(“)q2 

(o) = -. 

na=i ch(g - Pa^) ch{q + /ii“^) 

If it terminates with an antisymmetric hypermultiplet we have 


(2.45) 




327riM(“)p 

chp 




327riM(“)p 

chp 


^-TTQ 


(2.46) 


where is the antisymmetric mass. 

Note that the expression 02.431) assumes there is at least one U{2N) node. There are 
two relevant cases without U(2N) nodes: For the single node Sp{2N) theory with an 
antisymmetric hypermultiplet the density operator is given bj@ 


^ / ^27TiMp ^27TiMp 

p = - sh.(2q)F^^’ (q) ( sh q —^-ch g + ch q —^-sh q 

2 \ cap cap 


(2.47) 


For the Sp{2N) x Sp(2N) theory the density operator is given by 


g27rjm(®)p ^ ^iximF'lp 


P = 


chp 


■ sh(2g).F*^°^(g)- 


chp 


sh(2g).F*^^^(g) . 


(2.48) 


3 Mirror symmetry 

In [13] it was found that pairs of mirror dual A-quiver theories give rise to Fermi gas 
formulations which are related to each other by simple linear canonical transformations 
acting on their density operators. In this way, once one understands how to go from a 
quiver theory to the Fermi gas density operator, one can with little effort hnd the mirror 
map between dual theories. In this section we show that this continues to hold true 
for the mirror pairs involving iD-quivers and linear quivers. In this section we consider 
theories with mass and FI term deformations, since those get mapped to each other under 

^^This theory was previously studied from the Fermi gas perspective in m, where a rather different 
density operator was obtained. We compare the different formalisms in appendix iDl 
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mirror symmetry. For simplicity we restrict our attention in the remainder of the paper 
to theories without Chern-Simons terms 

We proceed with a number of examples, starting in each case with a Z)-quiver, and 
demonstrating how the canonical transformation 

q^-p, (3.1) 

maps the density operator onto that of the mirror linear quiver theory. 


3.1 Zl 4 -quiver with two fundamentals 

The hrst example we consider is the T) 4 -quiver with a fundamental hypermultiplet on 
two of the terminal U{N) nodes. This example is somewhat special as there are three 
inequivalent ways of pairing up the terminal U{N) nodes, which leads to different density 
operators. The canonical transformations of these three descriptions are related to three 
different mirror theories. The existence of several mirror dual theories was already noted 
in section 4.4.3 of m- 



Figure 6 : T) 4 -quiver and its mirror dual theory with Sp{2N) x Sp{2N) gauge group. 


The hrst possibility is to pair the two terminal nodes without fundamental matter 
and the two with fundamental matter, as shown in hgure [ 6 l Within each pairing we 
distinguish the FI parameters of the two U{N) nodes by giving one of them a prime. We 
do not turn on masses for the two U{N) fundamental hvpermultipletsF^ The density 
operator of the 1)4 theory can be read from 02.421) 


1 1 


P = 1 


4 chp 


g27riC(0)g^J)Tg27riC'(0)5 ^2iTiC'<-°'>q ^^P ^2TTiC<-°'>q 

chp chp J clip 


,2Tri(Wq 


1_/^^27riC(2)g^J)T^27riC'(2)g ^2TriC'<-^)q P ^27riC^^'lq\ _ ^^27ri(:W 


chp chq 


chp 


chp 


(3,2) 


chq chp 


To map this to the density operator of the mirror dual theory, we hrst use the relation 


- /(P - C). 


(3.3) 


to simplify the terms in parenthesis 

^ni{C-C)q^^-ni{C-C)q + g-«(C-C')9 

ch p ch p 

_ sh(p + IC - K) sh(P - + K) _ _ 2 sh 2 p _ 

ch(p + ic' - K) - K' + K) ch(p + K - K') ^^(p “ ic + |C') ’ 

^^One of these mass parameters can be removed by shifts of the matrix model eigenvalues. The other 
is mapped under mirror symmetry to a “hidden” FI parameter 113117], which does not have a clear 
interpretation in the mirror gauge theory. 
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By further commuting exponential factors we get 


p = ^ _ sh 2p _1_ 1 

ch {p + ch {p + ch {p — ch (p — ch {p + /ig^^) 

ch q ch (p + /2(°)) ch (p — p(°)) ch q ch (p — 

with 

^(1) _ 1(^(0) + (.,(0)) ^ ^(1) _ l«-(0) _ ^,(0), _ ^(1) _ _^(1) _ 1(^,0) ^ ^,(0)) _ 

yo)_ 1(^(2) _c-(2))^ (3.6) 

ft_-C<»-i(ci“> + C''“> + C'"> + C''"’). 


Now we can act on the density operator by canonical transformation fl3.ip . In addition 
we conjugate the operator to remove the exponential factors at the beginning and end, 
which clearly does not alter the spectrum. This gives 


^ 2 nifhp gp2g sh2g 

chp ch (g + ch (g — chp ch (g + fla^) ch (g — fia^) 


(3.7) 


As advertized we recover the density operator for a linear quiver with two Sp(2N) nodes. 
One with a single fundamental hypermultiplet and the other with three 02.481) . The 
relations between the FI and mass deformation parameters of the mirror dual theories are 
expressed in 03.6p . where m is the bifundamental hypermultiplet mass and are 

the fundamental hypermultiplet masses of the dual theory. This mirror map generalises 
slightly the one found already in [6], allowing for and ^ which 

translate to additional mass deformations in the dual linear quiver theory. 




Figure 7: l) 4 -quiver and its mirror dual theory with U{2N) gauge group and antisymmetric 
hypermultiplets. 


A second way of obtaining a density operator for the D 4 -quiver comes from pair¬ 
ing nodes with and without fundamental hypermultiplet, as in hgure [71 After shifting 
eigenvalues to remove masses we are left with only one nonzero mass /i for one of the 
fundamental hypermultiplets. From fl2.42p we can write down the density operator 


1 1 




P = 1 


4 chp 
1 




2'iTil^^^') c, 


^2TTi(Wq 


chg chp chp chg y chp 


(3.8) 


chp \ ch(g -I- fi) chp 


^2TTi(Wq 


chpch(g + /i) I chp 
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Once again we start by manipulating the expression of the density operator, using the 
relation fl3.3jl 


p = 


4 ch(p + ch(p - 


27ri(c('’)+c'<°4c<ibg sh(p + + —) sh(p - (4°^ - 


+ 


ch q ch{p + ^4°) + ch(p — 


ch q 




ch{p + 4^) ch(p - 

.2«(c(2)+C'(^)+C(i)). sh(p + C'(2) + sh(p - C'(2) - e2«(C(^)+C'(^)+C(i)).\ 

ch(g + /x) ch(p + C'(2) + ch(p-C'(2) - ^) ch{q + p) j 


X e 




We now use the identity 


^27ri(^q 


sh p + C' sh p — e 


2Tri(q 


+ 


ch(g + p) chp + (' chp — (' ch(g + p) 




^27Tifip 


^27ri(^q 




ch(p — (') ch g 

to bring the density operator into the form 

^ ^g-27ri(C(^)+C^^^+C'^^bMg-’r*C^^''Jg27ri/xp 

^ 4 

-2iTiiip 1 


g^p _ 


p2TxiC,q 

chg 


-Tip 


^-27vipp 

ch(p + (') 


ch (p + ch (p + ch (p — ch (p — ^4°) — 


,27ri(C(°)+C''°4c(lbQ 




ch q 


+ e 


— TVp 


j27ri(C(“)+C'™+C^^bQ 
ch q 


„-7rp 


^2wifj,p 


ch (p + ch (p + ch (p — ch (p — — C'(2) j 


^2ni{C(^)+C^^)+(^^))q 


o'^Pl 


ch q 


-e^P + e-^P- 


,2«(c(2)+c'(2)+C(i))q 
ch q 


g-7rp I g-27riMPg7riC(l)|J 


(3.9) 


(3.10) 


(3.11) 


Applying the canonical transformation fl3.ll) and removing the exponential factors at 
the beginning and end by conjugation, we obtain (up to an overall phasethe density 
operator of the (second) mirror dual theory, which is a U{2N) theory with two anti¬ 
symmetric hypermultiplets of masses Mi, M 2 , four fundamental hypermultiplets of masses 
Pa, CK = 1, • ■ ■ ,4 and an FI parameter The explicit mirror map between parameters is 
given by 

Ml = -((0) _ ^'(0) _ ^(1)^ = _^(2) _ ^/(2) _ ^(1)^ 

Pi = h = /^3 = /^4 = 

C = h- 

^^Such phases are unphysical, see footnote [TUI 
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Finally we may think of other ways to pair the U{N) nodes of the iD 4 -quiver which 
are similar to the two cases above. For instance we can consider the exchange 
in fl3.8p . This symmetry, which is completely trivial on the D-quiver side, manifests as a 
relation between two mirror U (2iV) theories which differ only by the values of their mass 
parameters 

Ml —> Ml + /ii + /i 4 , M2 M2 — fli — , 

ftl —> —fli, fl2 fj-2, P-S ^ h3) , (3.13) 

C-C- 


3.2 Ds-quiver 

The next example of mirror map involves a Hs-quiver with n fundamental hypermultiplets 
on one U{2N) node and a single fundamental hypermultiplet on one U{N) node. The 
quiver is shown in hgure [HI Shifting eigenvalues to set all bifundamental masses, and one 
of the U(2N) fundamental masses to zero, we hnd the density operator (with obvious 
notations for the mass parameters) 


p = U_ I g27riC(“)q^ihPg27riC'™Q + ^2-KiC^°)q 

4 ch p \ ch p ch p y ch p 


,27riC(l)i5 






27riC'(3)g 


clip ch q Yllzl ch {q + chp \^ch (g + chp 


(3.14) 


+ e 


2’jTiQ'^^^q 


,shp 




^2TTi((^\ 




chp ch (g + p^^^) ) chp ch g n”l} ch (g + pL^^) chp 




Figure 8: Zls-quiver and its mirror dual theory. 


Using the identities fl3.3p . fl3.4p and fl3.10p . we can write the density operator as 
2 chp ch(p + ch(p — 


X g^4C^°UC'^°bq ^ ^2nicmq 

chp 


^ g27riC(2)q ^ 

chp ch g Hall ch (g + pi^^) chp 


g27ri/j^®^p 

ch(p — C'1^1) 


2ni{(:(3)+C<^3))q 

gTi-p_gTrp 

ch g 


27ri(C(3UC'(®b'? 

+ ^- 

ch g 



g-27rip^®^p 

ch(p + 


X 


X e27riC(2)q ^ 

chp ch g Ylal\ ch (g + pL^^) chp 


g27riC(l)q 


(3.15) 
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Further commuting exponential terms using fl3.3jl we obtain 


p = ch [p + ~pf) g- 2 «(cw+c'( 3 ))^( 3 ) 


^—2iTirhq 


sh 2p 


^—27virhq 


n?=i ch (p + pf ) ch (p - pf) ch q 


'n-l 

'.a=l 


^27vi^^^'^p 


ch q 


,27riC(")j 




^—2'KiMq 


chpch(p + p"') \ chg 


g^P + g-^P 


^—2'KiMq 

ch q 


3-2vriC(")p 


-7rp 


ch p ch (p — ) 


' n—1 


nchg"" 

^a=l ^ 


j!_„-27riC("-“)p \ ^ _^ 


chg ch(p + p®) 


-7rj(C(°) + C'(°b9 


with 


(3.16) 


in - -pl»' - - C'^>. M - -e' - C'<=' 

(0, _ . ^,0, _ 1 _ -((,) _ . ^,0, ^ _ - (0, _ 1 ^ 1 ^ _ 


n 

^(») = _^'(3)_ 


((“) _ - pp - 0-2, 


(3.17) 


Applying the canonical transformation fl3.ip and removing terms at the beginning and 
end of 03.161) by conjugation, we obtain up to an overall phase (see footnotePIO]) the density 
operator of the mirror theory. This is a linear quiver with one Sp{2N) node with three 
fundamental hypermultiplets, connected to n U{2N) nodes, where the last U{2N) node 
has one antisymmetric and two fundamental hypermultiplets. By shifting eigenvalues the 
masses of hypermultiplets transforming in the bifundamental representation of two U{2N) 
nodes and the mass of one of the U(2N) fundamental hypermultiplets can be set to zero. 
This leaves us with the masses: fh for the Sp(2N) x U(2N) bifundamental hypermulti- 
plet, M for the antisymmetric hypermultiplet, ftfl 3 for the three Sp{2N) fundamental 

hypermultiplets, for one U{2N) fundamental hypermultiplet. Moreover the theory 
has FI parameters for a = 1, • • ■ , n. The mirror map between parameters is given by 
dSTH). 

This general approach to hnding mirror maps can be easily extended to more general 
D-quivers. For mirrors involving a ‘bad’ linear quiver, proposed in m the matrix model of 
the linear quiver is divergent 0 while the matrix model of the D quiver is still hnite. As we 
mention in the discussion section, the density operators of the ‘bad’ linear quivers can be 
matched to those of the ‘good’ D-quivers with a simple replacement, alas, this replacement 
is rather ad hoc and we do not know whether it represents a true regularization of the 
non trace class density operator of the ‘bad’ quiver. 


4 Computing the partition function 

Having re-expressed the partition function of D-quiver theories as a free fermion par¬ 
tition function, we proceed now with its evaluation following the technique developed 

was proposed in that the divergence of ‘bad’ quiver matrix model is related to a mismatch 
between the R-symmetry group in the UV localization computation, and the R-symmetry group at the 
infrared fixed point. 
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in [15]. More precisely we compute the perturbative part in the large p expansion of the 
grand potential J and extract from it the perturbative part in the large N expansion of 
the partition function Z, which turns out to be the Airy function fll.hp . Our strategy to 
compute J has new ingredients compared to [15], in particular we use a simplifying re¬ 
cursion method to compute the perturbative expansion of Tr p\ for arbitrary / ^ 1. This 
method can be used for any density operator p. We also rely on the integral representation 
of J, following [IH], to evaluate its perturbative part. Masses and FI parameters introduce 
extra difficulties in the computation and we set all of them to zero in this section. 

We proceed to present this strategy and provide the exact evaluations for all the 
D-quivers considered in this paper. 


4.1 General strategy 

The starting point for our computation of the perturbative part of the D-quiver (and 
linear mirror) partition functions fl2.25p is the Fermi gas reformulation 


Z(N) 



(4.1) 


where we recall that counts the number of cycles in a. The standard analysis [15] is 
to factor the integral into its closed loops 


J' dXi ■ ■ ■ dXi p(Ai, A 2 )p(A 2 , A 3 ) • ■ ■ p{Xi, Ai) 


Trp^^Zi. 


(4.2) 


These loops of course correspond to the cycles of the permutation cr, and so the summand 
of fl4.ll) depends only on Zi and the conjugacy class of a. Conjugacy classes of Sn can be 
labelled a set of integers {mi}, where mi is the number of cycles of length 1. In terms of 
this labelling we have 

2’T'o- i i 2™'i ' ' 

and the number of permutations in a given conjugacy class is given by 


A^! 

Yli 


(4.4) 


With these combinatorics fl4.ip becomes 


zw-Jl'Yl 

{mi} I 


i^Zi)^^ (_l)h-l)mi 
mil 


(4.5) 


where the primed sum denotes a sum over sets that satisfy n, Im, - N. 

The computation of Z{N) thus boils down to the evaluation of Zi. To compute results 
for large N, the standard approach [15] is consider the grand canonical partition function 

00 

5(2) = 1 + 2 ] Z{N)z^ = z = e\ (4.6) 

v=i 


where p is the chemical potential, and J(/i) the grand canonical potential, given by 


Jill) 


00 


-T 

1=1 


{—lyZie^'- 
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(4.7) 
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The strategy is to first find an expression for Z;, then to resum the expression and obtain 
J(/i) using fl4.7p . and hnally to recover Z{N) by computing 


Z(N) 


1 

27ri 


J ^Q—7Zi 




(4.8) 


where /io can be chosen arbitrarily without affecting the result. 

In practice, computing Zi exactly for arbitrary I is highly non trivial. To make the 
problem tractable it is useful to reformulate it within Wigner’s phase space [15]. For a 
general review of Wigner’s phase space see [50] ; here we simply summarise the properties 
that we require. 

The Wigner transform of an operato]0 A is given (with h = by 





(4.9) 


Some important identities are 


{Ab)w = Aw * Bw, 


* = exp 


i 

47r 


dqdp 


dqd 



In the language of phase space Zi becomes 


Tr(y4) = J dpdqAw ■ 

(4.10) 




J' dpdq pw * Pw ■ ■ ■ * Pw ■ 


(4.11) 


We generate an expansion for the integrand of 04.111) by performing a derivative expansion 
of the star products 04.101) . To this end, we introduce into the star product an expansion 
parameter e which will be set at the end to 10 


★ = 


exp 




dqdp 


= 1 + 


le 

47r 


dqdp dqdpj 


dldl + dldl-2dq,pdq,p 


(4.13) 


+ 0{eA 


In [T5], the role of the expansion parameter e was played by the Planck constant 
h, which was proportional to the Chern-Simons level of the ABJM theory. We do not 
have a tunable h here, however it still proves useful to consider the derivative expansion 
associated to e, as we now explain. 

The e expansion of the integrand of 04.111) takes the form 

{p^)w{p,q) = 0 Pi{n){p^q)- (4.14) 

n^O 

avoid confusion between phase space variables p,q, and the canonical position and momentum 
operators, we give all operators hats for the remainder of this section. 

^^In addition to the differential expression for the star product (Id.ldll . it also has an equivalent integral 
representation 

(/ * '?) = ^ J dq'dp'dq”dp"f{p + p',q + q')g{p + p”, q + q")e*^PP^'p"-A p') . (4.12) 

In some cases (for instance when f or g involve delta functions) (14.1211 produces an exact result that 
is non perturbative in e m- In these cases extra care should be taken, because the perturbative star 
product (14.101) is not valid. 
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Note that the e factors come from the expansion of the star products present in as 

well as those arising when replacing the density operator fl2.42jl by its Wigner transform. 
Zi can then be evaluated order by order in e 


£ ^l(n) ) ^l(n) I dpdiQ Pl(^n) ■ 

n^O d 

Resumming each term using 04.71) then generates an e expansion for J(/i) 

/ _IV 7 


(4.15) 


n^O 


1=1 


(4.16) 


As for the A-quivers [15], we anticipate that for D-quivers (and their linear mirrors), 
J(yu) also admits an asymptotic expansion of the form 

j(^) = ^^3^B(e)p + A(e) + 0(e-^^), a>0, (4.17) 

where each of the coefficients A, B and C are given by power expansions in e. In principle, 
to obtain a meaningful result for A, B and C we should now compute and then resum 
an inhnite series of corrections in powers of e, which are really all of the same order since 
e = 1. From the study of A-quivers, it is expected that the expansions of C{e) and B{e) 
truncate at orders and respectively, so that the hrst few orders of the e expansion are 
sufficient to compute them exactly. We give a proof of this truncation for l)-quivers with 
equal number of fundamental hypermultiplets on each pair of terminating U{N) nodes in 
appendix [El We assume that it holds for the other quivers as wellEl 

It remains to plug the result 04.171) into 04.81) to extract the perturbative part of Z 
at large N. In practice the evaluation is done by setting the contour integral parameter 
Po to the saddle point p* of the integrand and by extending the contour along all the 
imaginary axis to {p* — zoo, p* + zoo). As explained in [T^[22] . this change of contour does 
not affect the perturbative part of the result. The integration leads to the Airy function 
behaviour of the partition function at large N, which is our main result 


Z(N) 


C->e'^ Ai 


C-(N-B) +Z„p(Af), 


(4.18) 


where Z^p{N) denotes non-perturbative, exponentially suppressed contributions, and we 
note that the undetermined coefficient A only affects the overall prefactor. 


4.2 Recursive formula for {p^)w 

In this subsection we present a simple recursive approach for evaluating the coefficients 
in the e expansion of (p*)w 114.14p . This comes from the e expansion of 

= {p^)w * Pw ■ (4.19) 

One hrst needs to evaluate the e expansion of pw (which is due to replacing all operator 
products by star products) 

Pw (p, g) = 2 e" p(n) {p, q) , (4.20) 

n^O 

^®We checked that C and B do not receive contributions at order for all IJ-quivers. 
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which also serves as the initial conditions for the recursion pK^n) = P{n)- 
At order e°, equation fl4.19p gives 


which is solved by 
At order we get 


Pi+i{o) — Pl{0)P{0) 


Pi{o) ~ P{o) ■ 


Pi+i{i) ~ P{o)P{i) + P{o)Pi{i) ^ Pi{i) ~ ^P(o) P{i) ■ 

At order we then get 

_ , l I T 1—1 2 ^ l—2j c\ ( // •• •/2 'v 

Pi+i( 2 ) - P(o)Pi{ 2 ) + P{o)P{ 2 ) + ‘P(o) P(i) ^ 32^^(o) ''['^\Pio)P{o)P{o) ~ P{ 0 )P{ 0 )) 

+ - 1)(P(0)P(0) + P(0)P(0) ^ ‘^P{0)P{0)P{0)) 

where 

/(P, q) = dp f{p, q) , f'{p, q) = dq f{p, q) . 

Solving the recurrence relation, with initial condition pi( 2 ) = p( 2 ) yields 

Pi{2) = lp\o)P{2) + - 1)P(0)V(1) - ~ 3(P(0)P(0)P(0) - P(0)P(0)) 

+ (/ - 2)(p(0)P(0) + P(0)P(0) ^ ‘^P[o)P{o)P[o)) 


(4.21) 

(4.22) 

(4.23) 

(4.24) 

(4.25) 

(4.26) 


This procedure can be straightforwardly continued to higher order in e. Finally we plug 
the expansion coefficients into fl4.15p to obtain the e expansion of Zi. In particular this 
gives 

Zi{o) = J dpdqp[o) , 

Zi{i) = J dpdqlp\Q)P^^) , (4.27) 

Zi( 2 ) = J dpdq (^p\q) P{ 2 ) + “ l)P(ofP(i) “ “ ‘^)p[of Plo)Pio^ ■ 

Where the last line follows from integrating fl4.26p by parts, and using the fact that for 
A OT D quivers we always have the decomposition p^^-^ = t{p)u{q). 

We stress that this algorithm is very general and can be applied to any A or D-quiver. 
All one needs in order to compute Zi is to plug into fl4.27p the e expansion of the density 
operator itself. 

4.3 Computing Z(N) for D-quivers 

We now show how to apply the approach outlined in the previous sections and compute 
Z{N) for a generic D-quiver with arbitrary number of nodes and arbitrary number of 
fundamental hypermultiplets on each node. We set here all masses and FI parameters 
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to zero, as this simplifies the explicit evaluation of the phase space integrals. The quiver 
diagram is shown in figure [H and the density operator is given by fl2.42jl with 


a) = - 

ch q 




1 n'(“) 

cn q 


(4.28) 


We first work out the e-expansion of pw itself, which is then plugged in to the result 
of the recursive formula 04.2711 . 

Using manipulations similar to those used in section [3] 


shp 


-h 


shp 


ch"'°\ dip ch"''% ch”''% dip ch"'°\ 


^j^min(n(0),n'(0)) ^ chp 


★ ( sh p ★ 


^]^|n(0)-n'(0)| ^ 


★ ch p -I- ch p ★ 




★ shp 


1 


1 


Chp cJ^min(n(0),n'(0)) ^ 

1 


1 2 sh 2p 1 

★ - - m - TTTs. - 


(4.29) 


^j^min(n(0),n'(0)) ^ chp g chp (.J^min(n(0),n'(0)) ^ ’ 

where in the last step we evaluated the expression inside parentheses using the exact star 


product T :.12p . This allows us to write the Wigner transform of the density operator 

02:4211 a^pl 



= e'r(p) * * e^^Pd2Uoi,) ^ ^T{p) ^ 


A-l 

n, 

fc=i 


^ gT(p) 


★ e 


S(p)+2(7A(g) 


A-2 


* Yl 


fc=i 


(4.30) 


where 


and 


S{p) = logsh2p, T(p) = log 


chp’ 


Ui{q) = log 


ch’'* q ’ 


(4.31) 


Po = I 


Px = I 


Pi = min(?7,^°\ 


Pa- 1 = min(n*''^\ 


Pi = 7 = 2, • • • , L 


A = L + 2. 


(4.32) 


A first useful manipulation is to conjugat^ the density operator into a more sym¬ 
metric form that resembles a palindrome 


PlY % V e^^Pd2U\{g) ^ (Fiv) * 



^ ^S{p)+2Uo{q) 


★ 



★ (dd) ^ x/(lS{p)+2U\{q) ^ 


(4.33) 


^®n* defined by ordered star multiplication 

Since we are ultimately computing Zi = Trp' we can use cyclicity of the trace to make conjugations 
p V~^pV, which in the language of phase space becomes pw {V~^)w * Pw * Vw 
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where ^ is the star square root, and its expansion gives 

7eS(p)+2c/,(g) = gis+r/, (^2su'' + 2SU'^ + S^U'l) + 0{e^)^ . (4.34) 

The reason for making this conjugation to palindromic form is that the e expansion of 
fl4.33p is purely in even powers of e. The easiest way then to compute the e expansion is 
to build it up from the central piece. We hrst compute 

gT ^ gS+2C/o ^ gT _ g2T+5+2C/o . (4.35) 


By plugging this result in, we can then easily compute 


Ui . ^ ^ S + 2Uo ^ ^ 


) 


e^* e *6 '"★e M-e 


_ g 2 T + S '+ 2 J 7 o + 2 J 7 i 


- f{U'' + 217'") - -^[/"(2T + 5 + (2r + ^)2) + C>(e^) 


87r2 ^ 167r2 

Continuing this procedure we hnd the full expansion of fl4.33p up to O(e^) 


(4.36) 


= e 2 AT + 25 + 2 SLot /. 


,2 A 


2 A—1 / k 

6 ^ NTH / 


87r2 


167r2 


k=0 \J=0 '^i=o 


2 + 2kf +{S + 2kf) 

=1 

16,r2 + +f^7d + 2Ar) 1+0(6*) 


(4.37) 


fc =0 


fc=o i=o 


The coefficients of this expansion are p(o) and p( 2 ) fl4.20p . which serve as the seed for the 
recursion and can be plugged directly into fl4.27p . We hnd that Z;(o) is given by 


Zi{o) = J dpdq P(o) = J dpdq 






A 

k=0 


(4.38) 


The expression for Zp 2 ) is considerably more involved, but it can be simplihed by 
integrating by parts to remove all double derivatives. Integrating by parts the hrst term 
in Z;( 2 ) fl4.27p gives 


dpdqlp\Q)p^2) 

\ rdpdg/2g2HAT+S+SLo^.) 
47r2 J 


2r(s+Ar)|;|;q(|;c/'-(2!/' 

k=0j=0 kj=Q 


i=0 


+ i S + ‘^kf - 21{S + XT)) 2 U) 


(4.39) 


fc=l j=0 

S(S + Ar)(( - 1) 2 U', - isls + 2\f )u'>) Yj q 

/ i=0 


k=0 
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The second (non zero) term in Zi( 2 ) gives 

I - 1)(^ - 2)P(o)V(o)P(o) 

J ^ \j=o / 


(4.40) 


Combining these expressions, substituting fl4.3ip and simplifying leads hnally to 
Z,(„ - |3(2( - 1)AI. - (41^ - ly 

+ 2th^p(^(4l^ - 1)(A - 1)^2 - 3An(/(A - 2) + 1) - 6Si) (4.41) 

- th^p(^(4/2 - l)(A - 1)^2 - 3AV2 + 3Ai/(2/(A - 1) + 1) - 12(^2 - Si)) 

A = ?7o + hA, = 


where 


4>J 


A (4.42) 

S 2 = ((A - +f + A^) - 2r72z(A - i). 

i>j i=0 

In order to evaluate the integrals appearing in Z;(o) and Zp 2 ), we require only the 
identity 


dx 


th“x (1 + e*^“)r (^) r (I) 
“ 2^+%Fp±|±i) 


Re(a) > —1, Re{b) > 0 


(4.43) 


Integrating 04.381) and 04.411) using 04.43p . simplifying by shifting the arguments of Gamma 
functions and choosing to express the result in terms of L = A — 2 we hnd 


1 r(/ + A)r(^/)r(L/) 

m 22«c+C7r3/2 r(z// + |)r((L + i)i +1) ’ 

7 r(/+ |)r(zz/)r(L0 

«(2) 3.22Zd+C+3 ^ ’ Y{iyl + |)r((L + 1)/ + I) ’ 


(4.44) 


where 

F{1) = 1/2(1 + (L + 1)(L + 2)) - 3Au - 3(Si + S 2 ) 

+ /(n2(L + 2)(L + 3) - 6 Az/(L + 1) - Q{L + l)Si + 6 S 2 ) - 2fiz‘^L{L + 1). 

(4.45) 

Having computed the e expansion of Zi, we can now compute the e expansion for J{p) 
by resuming each term using 04.161) . As for the A quivers, this gives some complicated 
hypergeometric functions [l5l|39], from which we can then extract the large p asymptotic 
expansion 04.17p . Ultimately, we are interested only in the perturbative part of this 
expansion. A recent paper [48] pointed out a very elegant and simple way to extract this 
perturbative piece of J(/x), just by evaluating a single residue involving Z;. 

The hrst step is to write a Mellin-Barnes integral representation for the inhnite sum 

dHl 

/_IV 7 pU rc+iX) tj 

/(n)M = -g --J |)^r(!)rH)Z„„,e", (4.46) 
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where c can be chosen arbitrarily in (0,1), and should now be regarded as a function 
of /, analytically continued to the complex plane. To see how this integral representation 
reproduces the inhnite sum fl4.46p . note that for p < 0 we can close the contour of 
integration around the region with Re(/) > c. Since itself has no poles in this region, 
the only enclosed poles are the simple poles of r(—/) for I = n e N'*'. Using the fact that 

(- 1 )” 

Resr(-0 = ^^, (4.47) 

l=n Til 


we recover fl4.46p . 

Since we are interested in the asymptotic region jj, » 0, we close the contour of 
integration around the region Re(/) < c. In this region there can be poles due to both 
Zi(^n) and r(/). The residue at / = 0 turns out to be the only one giving a contribution 
that is not exponentially suppressed at large /i. Therefore we can immediately evaluate 

r(/)r(-/)Z,(„)e'^ + 0(e-“^), a > 0. (4.48) 


Evaluating in this way the perturbative contributions to J(p) from T(o)(/^) and J( 2 )(/i), 
we hnd the asymptotic expansion of the form (14.17p . with C and B coefficients 


1 — 3Aiy + + L(3 + T)(z/^ ~ 1) ~ 3(Si + S 2 ) 

12Ln 


(4.49) 


As explained before, the coefficients C and B do not receive contributions from higher 
order terms in the e expansion, so that fl4.49p provides the exact result. These coefficients 
characterize the full perturbative part of Z{N) as an Airy function fl4.18p . up to the 
overall coefficient A which is undetermined by our analysis. 


5 Discussion 


One of the main applications of our results concerns holography. We have found the 
complete large N perturbative result for the partition function of (good) £)-quivers as an 
Airy function (14.1811 . The large N expansion of the free energy starts with the terms 


-logZ(iV) 



B 


AtH ... . 


(5.1) 


The leading term has the famous behaviour and its coefficient depends only on the 
parameter C. This coefficient has a simple geometric interpretation when the theory 
admits an M-theory holographic dual of the form AdS^ x SE 7 , where SEj is a tri-Sasaki- 
Einstein manifold. Holography predicts the relation isa 

c - JLvo 1 (S'£ 7 ) . (5.2) 

7r° 


The i)-quivers studied in this paper (with vanishing masses and FI terms) can be en¬ 
gineered as the low-energy limit of a stack of 2N M2-branes at an orbifold singularity 
in C^/Z 2 y X [l2l[T3], where 7 j 2 u and are discrete subgroups of SU{ 2 ) of type 
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A and D respectively. The additional data describing the distribution of fundamental 
hypermultiplets in the quiver should be encoded in four-form fluxes on vanishing cycles 
at the orbifold singularity [12], however the precise dictionary is not known. The mirror- 
dual linear quivers have the same M-theory holographic backgrounds. These M-theory 
backgrounds can be reduced to type IIA and T-dualized to type IIB in two ways by ex¬ 
changing the M-theory and T-duality circles, leading to two IIB backgrounds in A-dual 
frames. 

The backreacted geometry in the large N limit takes the form AdS 4 x l{'L 2 u ^ ©l), 
where the two orbifolds act separately on the two ^ SU{2) inside S'^. Recalling 
that the volume of S'^ is and the order of the 'L 2 V and quotients are 2v and 4L 
respectively, we obtain the holographic prediction 


C 


6 

TT® 24zyL 


1 


(5.3) 


in perfect agreement with fl4.49p . It would be interesting to hnd a similar nice geometrical 
interpretation in holography for the coefficient of subleading term Bj^/C in fIS.ip but this 
coefficient is very subtle and not even fully understood in the case of ABJM [S3]- 

A localization calculation in AdS^ BDl suggests that the result would be universal for 
any conformal 3d theory with an M-theory dual and enough supersymmetry. The D- 
quivers (and their mirrors) studied in this paper fall into this class, and indeed we found 
the Airy function behavior, like for the A-quivers. Note that this calculation does not 
predict the values of the C and B coefficients, rather they serve as arbitrary parameters 
of the calculation. 

In our analysis of mirror symmetry, we left aside the D-quivers whose naive mirror 
duals are bad linear quivers m- These are D-quivers with numbers of fundamental 
hypermultiplets in a pair of terminating U{N) nodes differing by two or more. The 
bad linear quivers have a divergent matrix model leading to a non trace class density 
operator. One example considered in m for instance, is the naive duality between a D 4 - 
quiver with two fundamental hypermultiplets on an external U{N) node, and an Sp{2N)‘^ 
linear quiver with four fundamental hypermultiplets on one Sp{2N) node, and zero on 
the other. Manipulations as in fl4.29l) give the Wigner transformed density operator of 
the l) 4 -quiver as 

1 sh2p sh2p . 

,2 * ,2 * 16 • (^•^) 

ch p ch q ch p 

The Wigner transformed density operator of the naive mirror dual linear quiver is 


—— ★ sh2g ★ 
chp 


1 

- ★ 

chp 


sh 2 q 
ch® q 


(5.5) 


these operators would be related by the canonical transformation fl3.ll) (and conjugation 
by if we also included the replacement (with n = 1 and m = 0) 


1 sh2g 1 sh2g 

ch"" p ch™' p ch” p ch^”’*'™ p 


This is an ad hoc regularization of the partition function, transforming a non trace class 
density operator to one that is trace class, based on the assumption that the regularized 
partition function should satisfy the naive mirror symmetry. In fact, the replacement 
fl5.6p (with suitable n, m) easily regularizes and identihes a mirror dual for any bad linear 
























quiver with no fundamental hypermultiplets on one or both terminal nodes, provided the 
theory has in total at least four fundamental hypermultiplets. 

It would be interesting to understand if fl5.6p can be derived from a proper regulariza¬ 
tion of the divergent integrals in the matrix model. 

One can think of several extensions to our work. Having found a Fermi gas formalism 
for the A and D quivers, it is natural to ask whether a Fermi gas formalism for the 
quiver theories of type exist. To study these theories one needs to hnd suitable 

identities to put the corresponding matrix models in the form of a partition function for 
non-interacting particles in Id. 

Further generalizations involve necklace quivers with alternating Sp and SO groups. 
Likewise one can consider linear quivers with terminating SO nodes or symmetric hyper¬ 
multiplets (see ini)- One could also try to hnd generalizations of the D-quivers, replacing 
one of the D ends with some other gauge groups or matter helds. 
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A Matrix models for quiver theories 

We provide in this appendix the factors appearing in the D-quiver and linear quiver 
matrix models. Each factor is associated to an A/" = 4 multiplet of the theory and is a 
function of the eigenvalues of the matrix model. 

For U{N) nodes we have the ingredients 


Eigenvalues: 
Weyl group order: 
Vector multiplet: 

Fundamental 

hypermultiplet: 

Antisymmetric 

hypermultiplet: 


Aj i = I,-- - , V, 
|VF| = V!, 


^vec = ]^sh^(Ai -\j), 


i<j 



i<j 


1 

ch(Aj -h /i) ’ 

1 

ch(Aj -I- Xj + M) 


(A.l) 


The expression for U(2N) is the same with the obvious change N —> 2N, and where we 
use the convension of upper case /, J for the indices. 
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For Sp{2N) nodes the analogous expressions are 


Eigenvalues 
Weyl group order 
Vector multiplet 

Fundamental 

hypermultiplet: 

Antisymmetric 

hypermultiplet: 


\i i = ,N, 

|VF| = 2^m, 


Zvec — ]^sh^(Aj — Aj)sh^(Aj + \j) PJsh^(2Aj), 

i 

1 


i<j 


7 fund 

hyper 


n 


yasym _ rT 
^ hyper | 1 


ch(Aj + p) ch(Ai - /i) ’ 
1 


n 


ch(V + \j + M) ch(Ai + \j - M) ch(V 


Finally for the bifundamental hypermultiplets we have 


1 

Xj + M) 
(A.2) 


U{N)xU{My. 
Sp{2N) X U{2N) : 
Sp{2N) X Sp{2N) : 


N M . 

nn_-_ 

ch(Ai - Xj - m) ’ 

n_._i._, 

\ } ch(Ai -Xj-m) ch(Ai + Aj + m) ’ 

nn—^^^— 

\ j + ch(Aj + Xj + m) ch(Aj — Xj + m) 


(A.3) 


B Combinatorics of the permutations Rt-^Rt 

In this appendix we show how we can simplify the partition function fl2.2ip 

2 (-1)" f n \Rr->Rr(k)) . (B.l) 

TGS2N /cG/C(r) 

by studying more closely the composite permutation Rt~^Rt and the set /C(r) of N 
integers in 1, ■ ■ ■ ,2N such that R{K,{t)) = /C(r) and Rt~^Rt{1C{t)) = /C(r). 

Let us label the N integers in /C(r) by /ci, ■■■, /cat, such that the action of Rt~^Rt on 
/C(r) can be represented in terms of a permutation a-r ^ Sn as 

RT~^RT{ki) = . (B.2) 

We can immediately see that Rt~^Rt acts on elements R{ki) in the compliment of /C(r) 
by the inverse permutation 

RT-^RT{R{ki)) = R{RT-^RT)-\ki) = . (B.3) 

This property means that Rt^^Rt is composed of pairs of cycles that take the forn:!^ 

{kik2 • • • ki){R{ki)R{ki.i) • • • R{k,)). (B.4) 

For a given a e S'at we can easily hnd r e S 2 N such that a = a,- by for example taking for 
each /-cycle in a the 2/-cycle in r 


ikiR{k)ki_,R{ki_,) ■ ■ ■ R{k^)). 


(B.5) 


Indeed there is a freedom in choosing the set /C(t) by including the elements of either one of each 
pair of cycles. 
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This particular choice of r is useful because T{ki) = R{ki), and it is made up of only- 
cycles of even length, so we can easily compute 

(-1)" n =(-1)'^", (B. 6 ) 

kelCir) 


where we recall that 


s{k) 


0, k = !,■■■ ,N, 

1, k = N + 1,- ■ ■ ,2N , 


(B.7) 


and Ut- counts the number of cycles in r. The second equality in fIB.bp follows because 
each cycle in r gives rise to a cycle in ar, and in the third equality we recognised that the 
expression appearing is nothing but the signature of 

Of course, there are many possible choices of r giving rise to the same Rt~^Rt and 
we should check that all of them reproduce fIB.bp . From any particular r we can generate 
all T giving rise to the same Rt~^Rt by taking 


r 


TFT . 


TT ^Rn = R 


(B.8) 


To solve this condition, vr can be any permutation of the form tt = 7ri7r2, where tti is any 
combination of the two cycles appearing in R (2^ possibilities), and 7r2 acts as two copies 
of some Sn permutation (A^! possibilities) 


7r2(^) = ct'(z), n 2 {N + i) = N + i = 1, ■ ■ ■ , N, a' e S 


N 


(B,9) 


It is easy to check that any deformation flB.SD leaves the left-hand side of (IB.6P invariant, 
and so the right-hand side indeed holds for any r. With this simplification we can rewrite 
the partition function fIB.ip as 


Z(N)- 


2^-«NP Z-i 


N 


{-ir 


(B.IO) 


i=l 


The summand of fIB.lOD depends only on the conjugacy class of (Jr, determined by the 
number of cycles mi of length 1. This means we can convert the sum over S 2 N permutations 
to a sum over conjugacy classes of Sn, if we know the combinatorics of the map from r 
into the conjugacy class of ar- 

We have already counted how many r give rise to any given Rt~^Rt permutation 
(2^A^!), so we just need to compute how many distinct Rt~^Rt permutations are associ¬ 
ated with a (Tr of a given conjugacy class. We know that Rt~^Rt generates all possible 
permutations made up of pairs of cycles as in fIB.dp . and we should count how many of 
them have mi pairs of cycles of length 1. Suppose we £x how we assign to ki, - ■ ■ ,kN 
and R{ki), • • • , R{kN) the integers 1, • • ■ , 2N. With this restriction we are just left with 
counting the number of ways to distribute the ki among the cycles, which is the same as 
counting the number of S'at permutations in a given conjugacy class 


A^! 




(B.ll) 


We then have the added possibility of generating more distinct Rt ^Rt (without alter¬ 
ing the conjugacy class) by exchanging ki R{ki). Note however that simultaneously 
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exchanging all of the ki within a given cycle leaves Rt ^Rt invariant, so this generates 
just an additional factor of 


2N 
2'n.^r ’ 


(B.12) 


Putting all this together and relabelling Xk^ 


Aj, fIB.lOj) becomes 


Z{N) 


± V tDl 

TV! ^ 2”'" 

ctsSat 


N 


d'^Anp(Ai,A.(,)), 


i=l 


(B.13) 


where we have absorbed the factor (IB.Ill) to promote the sum over conjugacy classes 
to a sum over Sn permutations. 


C Degeneracy of the spectrum 

In this appendix we show that for vanishing masses and FI parameters the expressions 
for the partition functions found in the main text in terms of density operators fl2.25p 
(which is identical to (IB.13P above), correspond to the partition function of fermions on 
a semi-inhnite line fl2.27p . We do this by proving that the spectrum of p splits into odd 
and even states whose spectrum is identical]^ 

The hrst statement amounts to p commuting with the reflection operator R which acts 
on states by i? |A) = |— A). For vanishing masses and FI parameters the density operator 
p in fl2.42p is given by a sequence of even or odd functions of p or g, with precisely two 
odd functions. Since Rf{q) = f{—q)R, Rf{p) = f{—p)R, for any function /, we see that 
R commutes with p. 

To show that the spectrum of odd and even states is the same, we prove that Tr(p^i?) = 
0 for all 1. Let us focus hrst on the case / = 1. We notice that the density operator takes 
the specihc form 

P = B^d (P) ( 1 )pBd^ (p, q) p ^, (C.l) 

with p a sequence of even functions of p or g, pi is the Hermitean conjugate and 

+ a = 0,L. (C.2) 

Using the propertied 

R^ = R, f{qy = /(g), /(p)l = f{-p ), (C.3) 

we can derive the chain of equalities 

Tr(p^) = Tr ((p^)l) = Tr (^pl) = Tr [RpB^i\p, g) pi S® (p, g)) 

= - Tr (b® (p, g) p (p, g) pi ^) = - Tr (p^) , 

where we have used {B^^\p, g))l = —B^^\p, g), the cyclicity of the trace and commuted R 
and B^^\p, g), producing a minus sign. This yieds Tr(p.R) = 0. The argument generalizes 
easily / ^ 2. 

Apart for a single even zero-mode, which is non-normalizable. 

^^Here f acts as transposition, since we consider only real operators. 
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To derive fl2.25p . notice that the effect of the projection in fl2.27p . fl2.28p is to add 
a factor of 1/2 to every cycle in a given permutation a e Sn, 


In 




(C.5) 


The same results hold for density operators of linear quivers fl2.43p at vanishing masses 
and FI parameters. The arguments are the same except that one must consider the 
momentum space basis \p), with for instance Tr A = ^dp(^p\ A \p). 


D Relation with previous Fermi gas formulation of 
Sp{2N) quivers 

In this appendix we compare our results to those in [TT], which also developed a Fermi 
gas approach to Sp quivers. As we show below, the two formulations are equivalent, but 
in the final expression for the coefficients B in the Airy function, we disagree with their 
results and explain why. 

The overlap between theories discussed in jTT] and our work is the single node Sp 
quiver with an antisymmetric hyper mult iplet and n fundamental hypermultiplets. The 
partition function is given by the matrix model 


Z{N). 


r " 

k^n 

•d i=i 


sk 2\i ch2Aj ni<j sk(k — k) + k) 


qiViV! ch^- A, R,, ch(k - A,) ch(A. + A,) 

In |TT] the matrix model flD.ip was manipulated with a modified Cauchy identity 

Y\i<Axi - XjAvi - VjAxiXj - lAviVj - 1 ) 


(D.l) 


+ yj){xiyj + 1) 


T(-irn 


^Xi + ?/cr(2) ) (^22/cr(i) 1) 


This identity, with the usual the usual replacements Xi 




gives 


kv! 2 

aeSiv i=l 


Z(N) = 

1\ I 

ctgSn’ 

Using the relationl^ 


nisk2Ai ch2Ai 


Yl^ ck” Xi ch(Ai - A^(i)) ch(Ai + A^p)) 


(D.2) 

(D.3) 


ch vrAi chA 2 . .l + i?, . 

ch(Ai — A 2 ) ch(Ai + A 2 ) ^ 2chp ^ ’ 

sh Ai sh A 2 _ / \ k “ ^ I \ \ 

ch(Ai — A 2 ) ch(Ai + A 2 ) ^ 2chp ^ 

from which one immediately has the corollary 

1 — R shql + Rshq 

2 ch p ch g 2 ch p ch g ’ 


(D.4) 


(D.5) 
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Note that R commutes with and so we can write ) which would otherwise be 

cn p ’ 2 cn r) ' chv 2 ^ 


ill defined 


2 chp 


chp 2 
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they rewrite the partition function in two equivalent forms 


Z(N) = 


I , K « 

(TsSiv *=1 *=1 


P± 


1 + -R 




(D.6) 


where 


sh^ g ch 2g 1 
clR”g chp ’ 


P- = 


ch 2q 


(D.7) 


ch gchp 


Indeed, this looks very similar to our rewriting of the matrix model as fl2.27p or fl2.28p 

N N 


P\ ^ I IV(.)> 


creSjv i=l i=l 

but with a different density operator 02.471) 


lsh2g/ 1 1 , 


(D.8) 


(D.9) 


All those expressions are in fact equivalent, since we can show that the spectrum of 
even states of p agrees with that of and the spectrum of odd states to that of p_. 
Indeed the projected operators are similar to each other by the simple manipulations 


1 — R 1 sh 2g / 1 

= 0-7^1 shg 


2 ch ""q 
1 sh2g 


2 chp 


R , 1, ^ R u 

ch g + ch g—-— sh g 


2 chp 


2ch2"g V'^“''2chp 
ch 2g ch g 1 — R 


, 1 — R , ch^ g 1 — R 

sh g—-— ch g H- - --— ch g 


shg 2chp 

R 


(D.IO) 


ch g 2 ch p 


chg = -^( p_- 
ch g \ 


ch g 


where in the hrst line we used that R commutes with even and anticommutes with odd 
functions of g, and in the second line we used OD.Sp . Similar manipulations lead to the 
relations 


1 

ch g 



chg 


1 ( l + R \ 
sh g 2 ) 


sh g 


chg / 1 + R\ 
ch 2g sh g y 2 ) 


ch 2g sh g 
ch g 


(D.ll) 


We can also derive the Airy function expression for these theories based on our analysis 
in sections E] and 01 The mirrors of this class of theories are D quivers with (n — 1) 
U{2N) nodes and and a single fundamental hypermultiplet on one U{N) node. The 
asymptotic expansion of the grand potential fl4.49p for this theory (L = n — 1, z/ = A = |, 
El = S 2 = 0) is given by 


C 


1 

27r2(n — 1) ’ 


B 



(D.12) 


The coefficient C is the same as found in [H], but B is not. The reason for the 
discrepancy is that in our formulation, the operator p has degenerate odd/even spectrum. 
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as shown in appendix O While the odd and even spectra agree with p+ and p_ of BB, 
those latter operators do not have degenerate spectra. The saddle point calculation of HU 
based on the Fermi surface of p_ considered the full spectrum of this operator and divided 
the end result by 2. Since the operator does not have a degenerate spectrum, this is merely 
an approximation that is good enough to evaluate the leading order term C, but fails for 
the subleading coefhcient B. 


E Truncation of the e expansion 


Here we show that the e expansions of C and B coefficients in fl4.17p truncates at order 
e° and respectively, adapting a similar proof from [15] for A quiver theories. 

We recall that corrections to (7, B and A coefficients at order can be computed 
from a single residue involving 

J(/x) = 2eV(„)(p), J(„)(p) = -^Rj3sr(/)r(-/)Zp„)e'^ + C>(e-“^), a>0. 

n^O 

(E.l) 

From this expression it is clear that if vanishes at / = 0, then can only correct 

the A coefficient. We should prove then that has at least a simple zero at / = 0 for 
all n > 2. 

It is useful to consider the Wigner-Kirkwood expansion of Zi |Tn|. The idea is to 
express Zi in terms of the Fermi gas Hamiltonian, Hw = — log* pw 

= S a = {{H - HwY)^. (E.2) 

r^O 'I 


We recall that the D density operator is given by fl4.33p 

= Y/e^^P)+‘iUx{q) * f n * ^s{p)+2Uo{q) 


A-l 


k=l 


(E.3) 


I I g^(p) ^ gUkiq) j ^ * Y (^S{p)+2U\{q) ^ 


Since flE.3p has an expansion in purely even powers of e, Hw likewise takes the form 

Hw = 7/(0) + e^77(2) + c^77(4) + • • • • (E-4) 

This gives 


r^O ■ 


dpdq e 


-IH, 


(0) 






1 - ; E E + 0(e) 9r . 




(E,6) 


We know that all of the integrals that appear in this expansion are of the form fl4.43p 

th“a; (1 + e'™)F (^) F (I) 


dx 


ch‘i 2''+iirr (!!±|±i) 


(E.6) 


The coefficients a and b are given by linear functions of I and so the integrals over p and 
q in flE.Sp can each produce at most a simple pole at / = 0, from the F (|). Therefore, 
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we need only concern ourselves with terms in flE.5jl with prefactors of at most order 
Discarding also (most of) the terms of order or less, we are left within 



For simplicity let us now restrict the density operator flE.Sp to cases with U^^q) = U\{q) = 
00 With this restriction, the exponentials in flE.Op can be freely exchanged with star 
exponentials, and we can straightforwardly evaluate the star logarithm by the star product 
version of the Baker-Campbell-Hausdorff expansion. The leading term gives 


e-^^(o) 


/ th^p V 

V ch^^ p ch^^ q ) 


(E.8) 


All e corrections to Hw are then given by nested star commutators involving T{p), S{p) 
and Uk{q)- An example of such a term contributing to if( 4 ) would be 


[T(p), [T(p), lT(p), [T(p), (7t(«)]]]]. = + 0(e«), 


(E,9) 


where we have used 


[f j 9 ]* = f * 9 ~ 9 * f = 9 } + higher derivative terms . (E.IO) 

ZTT 


As illustrated in flE.Op . in order to have a term at order e"' with only single derivatives 
acting on functions of p (or q), such a term has a single function of q (or p) with an 
derivative. 

This is important, because every term in the epsilon expansion has therefore at least 
one multiple derivative of iS, T or f/fc. From fl4.3ip it follows that these derivatives take 
the form 


A'A) - 


ch^ q 


th“g 

^ab 


asZ b^O 


ch^ q 


n ^ 2, 


(E.ll) 


and similarly for S, T with q ^ p. 

We now should combine these derivative terms with flE.SD in flE.7p . and integrate 
with flE.bp . It is clear that since the derivative terms contribute always a or the 
resulting Gamma functions can contribute at most a simple pole at / = 0. This guarantees 
the terms in (IE.71) with an P outside have at least a simple zero, leaving us with 


-/ 





(E,12) 


By the same reasoning, terms with multiple derivatives on both a function of p and a 
function of q have an overall which kills both of the poles one could get from 

^®The final term {Hw * Hw — H^) has still some piece, but as we shall see this also has at least a 
simple zero at Z = 0 

^®This corresponds to restricting the D quiver theories shown in figure[T]to a subclass with 
and nT) = n'T). We expect that theories outside this subclass also have B and C coefficients truncating 
at e^. We have verified that this holds true up to e^. 
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integrating. The remaining terms in flE.12jl which don’t obviously have a simple zero are 
those where all derivatives on functions of p (or q) are hrst order, like flE.9jl . But such 
terms can be integrated by parts, for instance flE.Qp would give 


, (E.13) 

Integrating by parts pulls down an additional factor of I, which guarantees that there is 
an overall simple zero at / = 0, since the integral on the right hand side still produces just 
a simple pole. Since we have shown that at order and higher Zi has at least a simple 
zero at / = 0, this concludes the proof that C and B do not get contributions beyond 
order e^. 
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